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ABSTRACT 

This Letter investigates the physics that is responsible for creating the current 
system that supports the outgoing Poynting flux emanating from the ergosphere 
of a rotating black hole in the limit that the magnetic energy density greatly ex- 
ceeds the plasma rest mass density (magnetically dominated limit). The under- 
lying physics is derived from published three-dimensional simulations that obey 
the general relativistic equations of perfect magnetohydrodynamics (MHD). It 
is found that the majority of the Poynting flux emitted from the magnetically 
dominated regions of the ergosphere has a source associated with inertial effects 
q outside of the event horizon. 
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1. Introduction 



There are two known theoretical mechanism for producing field-aligned outgoing poloidal 
Poynting flux, S p , at the expense of the rotational energy of a black hole in a magne- 
tosphere that is magnetically dominated. There are electrodynamic processes collectively 
called Blandford-Znajek mechanisms in which currents flow virtually parallel to the proper 
magnetic field direction (force-free currents) throughout the magnetically dominated zone all 
the way to the event horizon (Blandford and Znajek 1977; Phinney 1983; Thorne et al 1986). 
Therefore, these electrodynamic currents have no source within the magnetically dominated 
black hole magnetosphere. Alternatively, there is the GHM (gravitohydromagnetic) dynamo 
in which large relativistic inertia is imparted to the tenuous plasma by black hole gravity that 
in turn creates a region of strong cross- field currents (inertial currents), J x , that provide the 
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source of the field-aligned poloidal currents, J p , that support S p in an essentially force-free 
outgoing wind, i.e., V • J w d.J ± /dX ± + J p /dX p m at the source (Punsly 2001). The 
force-free electrodynamic current flow is defined in terms of the Faraday tensor and the cur- 
rent density as F^ V J U = 0. As a consequence, J ■ E = 0, so S p in a force-free magnetosphere 
must be injected from a boundary surface. The two types of sources associated with these 
two types of Poynting fluxes are quite distinct in the ergopshere: the inertial 

current provides the J • E source in Poynting's Theorem and the force-free (electrody- 
namic) component of S p emerges from a boundary source at the event horizon. 

Numerical models can be useful tools for understanding the source of S p emerging from 
the ergosphere of a black hole magnetosphere. Some recent three-dimensional simulations 
in De Villiers et al (2003); Hirose et al (2004); De Villiers et al (2005a,b); Krolik et al 
(2005) show S p emanating from magnetically dominated funnels inside of the vortices of 
thick accretion flows. In these simulations, J 1 - ~ J e (in Boyer-Lindquist coordinates which 
are used throughout the following) inside of the ergosphere, since the poloidal field settles to 
a nearly radial configuration early on in the simulation (Hirose et al 2004). In principle, one 
can clearly distinguish the amount of S p emerging form the ergosphere in a simulation that 
is of electrodynamic origin (as proposed in Blandford and Znajek (1977)) from the amount 
due inertial effects (the GHM theory of Punsly (2001)) by quantifying the relative strengths 
of S p emerging from the inner boundary (the asymptotic space-time near the event horizon) 
with the amount created by sources within the ergosphere. In the high spin rate simulations 
in questions, over 70% of S p emerging from the ergopsheric funnel is created outside of the 
inner boundary. 



2. Physical Quantities in Boyer-Lindquist Coordinates 

The Kerr metric (that of a rotating uncharged black hole), g^ u , in Boyer-Lindquist 
coordinates (r,9,4>,t), is given by a line element that is parameterized by the black hole 
mass, M, and the angular momentum per unit mass, a, in geometrized units (Thorne et al 
1986). We use the standard definitions, p 2 = r 2 + a 2 cos 2 9 and A = r 2 — 2Mr + a 2 , where 
A = at the event horizon, r + = M + V M 2 — a 2 . The " active" region of space-time is the 
ergosphere, where black hole energy can be extracted, r + < r < r s = M + \J M 2 — a 2 cos 2 9 
(Penrose 1969). 



- 3- 



2.1. The Toroidal Magnetic Field Density and the Cross-Field EMF 

The flux of electromagnetic angular momentum along the poloidal magnetic field direc- 
tion is the component of the stress-energy tensor, = [l/(47r)]i 7 a i 71ar , in the approxima- 
tion that the field is radial. In steady state, the electromagnetic angular momentum flux per 
unit poloidal magnetic flux is the toroidal magnetic field density: — B T = y/^ g F 9r , where 
g = — p 4 sin 2 # (Phinney 1983; Punsly 2001). Similarly, the electromagnetic energy flux 
along the poloidal magnetic field direction, S p , is the component, T t r = [l/(4ir)]F ta F ar , 
in the approximation that the poloidal field is radial (Thorne et al 1986). In steady state, 
the energy flux per unit poloidal magnetic flux is — (Q F /c)B T , where Q F is the field line 
angular velocity (Phinney 1983; Punsly 2001). Consequently, B T is useful for quantifying 
the energy and angular momentum fluxes as steady state is approached. From Ampere's 
law, 

V^J 6 = B T r + {y/=gF*)s . (2-1) 

At late times, as an approximate steady state is reached, one expects y^gJ e ~ B T r . There- 
fore, at late times, J e is a potential source for the current system that supports S p '. 

Even when a system has not reached a time stationary state, one can introduce a well- 
defined notion of Q F that becomes the field line angular velocity in the steady state. If 
the field is nearly radial one can simply define the expression, F te = —il F F d( p. With this 
definition, Q F is a function of space and time and in steady state it becomes a constant 
along a perfect MHD flux tube. Therefore, the EMF across the magnetic field is —£l F F e( p by 
definition. 



2.2. The Source of Poloidal Poynting Flux 

In Boyer-Lindquist coordinates, the curved space-time equivalent of the " J-E" source of 
S p is the term F ta J a . This notion is described by the integral version of Poynting's Theorem 
which is the integral of T" . v = F tu J u combined with Stokes' theorem (Thorne et al 1986; 
Phinney 1983) 

J \F t ^ + F tr J r + F te J e ] dV-jJ T/ dV = [1/(4tt)] j V=g~F ta F an dA; , (2-2) 

where dV = y^gdrd(f)d9 is the spatial volume element of a section of a thick spherical shell 
and — T t * is the energy density of the field and n is the normal direction to the surface area 
element, dA, of the Gaussian pillbox. In steady state, the source of S p is —il F F e<t) J 9 in the 
approximation of a radial field. Without the radial approximation, one needs to define a 
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poloidal field direction, B p , and a cross-field poloidal EMF, E± = —VL F B P then the source of 
S* p is E^J- 1 . This can all be setup with great mathematically complexity (see Punsly (2001) 
for sample calculations) and no additional physical insight. The reader should remember that 
the poloidal field is not exactly radial in the following and the simple notion of J 9 is used to 
approximate J L . 



3. The Source of Ergospheric Poynting Flux in the KDE Simulation 



The simulation "KDE" of De Villiers et al (2003); Hirose et al (2004); De Villiers 
et al (2005a,b); Krolik et al (2005) is of the most interest since it generates an order of 
magnitude more S p than any of the other simulations (Krolik et al 2005). The magnetically 
dominated funnel spans the latitudes 0° < 9 < 55° at the inner calculational boundary, near 
r + . Quantifying the magnetic dominance in the funnel is the pure Alfven speed, Ua = 
B p I (y/Annfic) , where n is the proper number density, /i is the specific enthalpy of the 
plasma and B p is the poloidal field strength. Within the funnel 10 < U\ < 10 4 . The funnel 
S p averaged over time and azimuth in KDE near r + is shown in figure 1. 

This contour map indicates a region of strong outgoing S p in the evacuated funnel, 30° < 
9 < 55°, r > r + . It is clear that S p suddenly diminishes close to r + at r ps 1.3M — 1.5M. 
Inspection of the contour map indicates that over 72% of S p is created within a thin layer 
near r 1.4M. Because of the saturation of the dark red color in the plotting routine, S p 
might be even larger above the switch-off layer than indicated in the contour map. Thus, 
we only have a lower bound on the strength of S p above the switch-off layer and it is likely 
that more than 72% of the energy flux is created within this thin layer. This effect even 
continues into the weak S p region closer to the pole at 20° < 9 < 30°, r > r + . 

In order to investigate possible source terms for S p , in (2.2), one needs to be explicit 
about what is plotted in figure 1. Whenever a time lapse of 80 M occurs within the high time 
resolution simulation, data is stored for a time snapshot. There are 75 of these snapshots 
that are averaged in figure 1, steps 26 through 100. Thus, the discrete time average of (2.2) 
is relevant , 
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(3-1) 



where the four sides of the Gaussian pillbox of integration in figure 1 are the four curves 
labelled "1 -4". The result of significance from this plot is that ^°° 26 J 3 y^gF ta F ar d9 d<p > 
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2 I Zl!=26 Ii \f-~9FtaF ar d6 d<p |, i.e. the time averaged S p increases across this thin volume. 
There are three possible source terms for this increase in S p when (3.1) is applied to the 
pillbox in figure 1. The time averaged field decay from a finite spatial region (the second 
term on the LHS) is not a physically viable source of long term S p (unless the simulation 
is pathological and keeps creating local field energy from numerical artifacts). The only 
reasonable choices are the J • E term and the conversion of latitudinal Poynting flux, S 1 , 
from the surface terms J^l^e J2+4 V~dFtaF ad dr d(p. We don't have time averaged plots of 
S ± , so we don't know how the magnitude of this term. Because of the symmetry at the 
pole, the only viable scenario is that S 1 - radiates poleward from its source in the funnel wall 
and somehow converts to S p in the pillbox. One thing that is unequivocal from figure 1, is 
that during the lifetime of the simulation, less than 30% of S p that reaches emerges from 
the ergospheric funnel came from the inner boundary, the rest was created external to the 
boundary. 

The existence of this putative source is a highly significant result and it is desirable to 
check it by other means. For example, since this occurs near the boundary how do we know 
that this is not an artifact of an inherent error in the plotting routine? Figure 7 of Krolik 
et al (2005) is a similarly averaged (in time and azimuth) plot of the magnetic angular 
momentum flux, T£ = [1 / (47i)]F^ a F ar . Figure 7 shows that T£ is created predominantly 
in the ergosphere, over 70% of the electromagnetic angular momentum flux in the funnel 
is created by sources within the ergopshere for 25° < 9 < 55°. The remaining fraction 
can be associated with electrodynamic sources (i.e., sourceless and emerging from the inner 
boundary). A consistent picture emerges from the time average of Ampere's law in (2.1) if 
the dominant source term for S p in (2.2) is F ta J a F te J e . 



4. The Cross-field Current Density in the Ergosphere 

In this section, we use a plot of the strong electromagnetic forces within the evacuated 
funnel to try to understand the physical mechanism that drives the source of S p in the 
simulation, KDE. Presently, there is no existing data that has been extracted from the KDE 
simulation that directly illustrates the electromagnetic forces or currents in the ergosphere. 
Fortunately, plots of the electromagnetic force already exist for the KDP model of De Villiers 
et al (2005a) which is characterized by a/M = 0.9. So, the best we can do is to explore 
electromagnetic force plots (such as figure 2) from a closely related simulation, KDP. This 
model has the second highest S p luminosity within the family of simulations, but it is an 
order of magnitude weaker than the KDE model. This should still be qualitative adequate 
for the following analysis, since it is stated in De Villiers et al (2005a) that the azimuthal 
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force exists in all the models in the same relative location, but it is strength correlates with 
spin and is the most pronounced in the KDE model. 

The evacuated funnel of KDP is dominated by magnetic energy, 10 4 > U\ > 10 (accord- 
ing to figure 3 of Hirose et al (2004)) in the ergosphere if 0° < 9 < 65°. The nonforce-free 
nature of the current density will be studied by means of a plot of the electromagnetic force 
(that were presented in De Villiers et al (2005a) and does not appear in De Villiers et al 
(2005b)), F*"J V = F&Jt + F^Jr + F^Jg, within the funnel at t = 8080M, averaged over 
azimuthal angle. 

One can scale the individual components of the strong azimuthal force that appears 
at r 1.5M in figure 2 if turbulence does not dominate the dynamics: F^Jt ~ a~ 2 F^ t , 
F^ r J r ~ a~ 2 F rt and F^ 9 J e ~ a~ 2 F et , where the quantity, a = V~K sin 9/^/g^, is the 
lapse function that represents a global redshift factor (Thorne et al 1986; Punsly 2001). 
The lapse function vanishes at the horizon and is therefore useful in expansions in a small 
dimensionless parameter near the horizon (at r=1.5M, or 2 = 58.8, 9 = 45°). These scalings 
yield the approximation 

2 Mr a 

F*»J V « F^J e « (fl, - fl) F^J 9 , (4-1) 

where VL = —g^/g^ is known as the ZAMO angular velocity and it approaches the horizon 
angular velocity deep in the ergosphere (see Thorne et al (1986) for more details). Figure 2 
indicates large changes in F^ v J v occur on the order of the grid size at r « 1.5M, ~ 0.01M. 
Thus, this must be a dynamic effect and does not derive from Va~ 2 . A consistent picture 
can be constructed by a large J e , enhanced by a factor of more than 20 compared to the 
upstream flow, that is responsible for the large force that initiates at a ~ 0.10 — 0.25. From 
Ampere's law in (2.1), this should be a prestigious source of B T and by (2.2) this should 
also be the strong source of S p , F te J e consistent with the analysis of the KDE simulation 
in section 2. Note that F$ a J a F^J 9 is the force that torques the plasma (Thorne et al 
1986; Punsly 2001). So, the Poynting flux generation in this scenario is associated with a 
strong electromagnetic torque. 

A physical explanation of the strong electromagnetic torque is the most basic concept 
of GHM. 

Namely, the field is rotating a rate that is slower than the enforced rotation velocity of 
the plasma that is dictated by the dragging of inertial frames in the ergosphere. Gravity 
tries to pull the plasma forward relative to the field as the horizon is approached, the back 
reaction of the field is to torque the plasma backwards back onto Larmor helices that thread 
the field lines. In the process, the field is partially overwhelmed by gravity and twisted 
forward azimuthally creating B T < even though the plasma is tenuous. Note that J 6 
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directed equator-ward in the northern hemisphere of the KDE simulation is of the correct 
orientation to both source B T < and to torque the plasma if B p > (for more details of 
this process see Punsly (2001); Semenov et al (2004)). The strong force occurs in figure 2 
throughout the range 10 < U\ < 10 4 indicating that the current is not driven across B p by 
the rest mass inertia, but by a relativistic effect imposed on the tenuous plasma by the black 
hole. Finally, GHM naturally explains the current closure within the global wind. Figure 
3 shows a remarkable agreement between the current system of the analytical GMH model 
depicted in figure 9.12 of Punsly (2001) and the KDE simulation. In the dynamo source, 
d(^/—gJ e )/d9 + d(y^gJ r )/dr 0, where J r is the field aligned current (dashed contours 
in figure 3) that support S p . 



5. Conclusion 

In this paper, we studied simulations of a magnetically dominated funnel of a rapidly 
rotating black hole, a/M = 0.998. A source that is responsible for creating over 70% of 
S p transported through funnel during the life of the simulation was found. Similarly, one 
can conclude that < 30% of S p emerging from the ergospheric funnel is from an inner 
boundary source, near the horizon. The small residual S p injected from the boundary into 
the accretion wind can be considered of electrodynamic origin. The distribution of S p in 
figure 1 is in contrast to the Blandford-Znajek solution in which essentially all the S p is of 
electrodynamic origin, i.e., it emanates from the horizon and passes through the accretion 
wind with minimal interaction, thereby maintaining a virtually constant value along each 
poloidal flux tube throughout the ergosphere. 

There were two possible sources for S p in the funnel, there are the GHM inertial cur- 
rents and S 1 - injected from the funnel wall. The latter would be a new source of S p in a 
magnetically dominated magnetosphere that has not been considered in previous literature. 
The putative S 1 is created by inertial effects in the funnel wall. We were unable to dis- 
tinguish between these two inertial sources with the available data. It was demonstrated 
that a GHM current explains the ergospheric sources of S p , the electromagnetic angular 
momentum flux, the large azimuthal forces seen in the ergosphere and global current closure 
in the wind zone. 

Further analysis of three-dimensional simulations about rapidly rotating black holes 
(a/M = 0.998) are needed to clarify the physics that creates S p . At least three consecutive 
time snapshots are needed in order to find the current distribution from Ampere's Law at 
each coarse time step data dump. Since the flow is highly turbulent near the event horizon, 
it would be important to then time average the current distribution. Higher resolution might 
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be required to resolve the issue of whether the source in figure 1 is from J E±, S or some 
unexpected numerical error. 

6. Figure Captions 

Figure I. The azimuthally averaged and time averaged (over 75% of the simulation that 
ends at t = 8080M) Poynting flux from the model KDE (a/M = 0.998) of Krolik et al (2005). 
The figure is a magnification of the inner region of figure 10 of Krolik et al (2005). It is an 
excision of a region, 0° < 6 < 65°, r > r + that is a little larger than the ergospheric portion 
of the magnetically dominated funnel, 0° < 6 < 55°, r > r + . The majority of S p switches-off 
in a thin layer near r = 1.3M - r = 1.5M (see Krolik et al (2005) for a description of the 
units on the color bar). This region is a source for the majority of the outgoing S p emerging 
from the funnel. A Gaussian pillbox is drawn as a dashed contour for use in Poynting's 
Theorem. There are 26 grid zones between the inner boundary, r = 1.175M and r = 1.5M. 
The plot is provided courtesy of John Hawley 

Figure 2. The azimuthally averaged azimuthal electromagnetic force in the evacuated 
funnel of the model KDP from the central region of figure 7 from De Villiers et al (2005a) at 
time t = 8080M in geometrized units. Magnitude information is expressed by color in code 
units (see scale above). The funnel wall is marked by the solid white line. The azimuthal gas 
pressure forces are much smaller in this region (De Villiers et al 2005a). Note that there is a 
strong electromagnetic torque between the boundary at r = 1.45M and r = 1.75M, an order 
of magnitude stronger than anywhere else in the ergosphere and the adjacent space-time. 
The plot only covers the funnel in the ergopshere in the restricted span, 45° < 9 < 65°, yet 
the full range of Alfven speeds is captured, 10 4 > U\> 10. The plot is provided courtesy of 
John Hawley 

Figure 3. An overlay of the time stationary current system from a GHM model repro- 
duced from figure 9.12 of Punsly (2001) and the Poynting flux distribution of KDE. The 
relative scales are set by the black hole (gray) radius, r + = 1.06M, and the inner boundary 
of KDE at 1.175 M. The GHM model assumes very different initial conditions than KDE, 
but the basic current topology should be common to all GHM magnetospheres: an enhanced 
strong cross-field current density (which increases equator-ward) in the ergosphere that is 
the source of a field aligned current system in a nearly force-free wind zone and a strong 
return current flow at the edge of the wind zone 
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